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Abstract 

A grand rotation curve of the Milky Way Galaxy is constructed, which covers a wide range of radius from 
the Galactic Center to ~ 1 Mpc, and is deconvolved into bulge, disk and halo components by least-squares 
fitting. We determined the scale radii and masses of the bulge and disk to be Mb = (1.652±0.083) x 10 10 M Q , 
a h = 0.522 ±0.037 kpc, M d = (3.41 ± 0.41) x 10 10 M Q and a d = 3.19 ±0.35 kpc. The dark halo was fitted by 
the Navaro-Frenk- White (NFW) density profile, p = p /[(R/h)(l + R/h) 2 }, and the fit yielded h= 12.5 ±0.9 
kpc and po = (1.06 ±0.14) x 1O _2 M pc -3 . The local dark matter density near the Sun at Ro = 8 kpc 
is estimated to be pf = (6.12 ±0.80) x 10~ 3 M Q pc~ 3 = 0.235 ±0.030 GeV cm" 3 . The total mass inside 
the gravitational boundary of the Galaxy at R ~ 385 kpc, a half distance to M31, is estimated to be 
M h ~ (7.03 ±1.01) x 10 n Af Q . This leads to the stellar baryon fraction of (M b +M d )/(M h + M d + M h ) = 
0.072 ±0.018. Considering expected baryon fraction in the Local Group, we suggest that baryons in the 
form of hot gas are filling the dark halo with temperature of ~ 10 6 K and emission measure ~ 10 _5 pc cm 2 . 
Such hot halo gas may share a small fraction of the observed X-ray background emission. 

Key words: galaxies: dark matter — galaxies: structure — galaxies: The Galaxy — galaxies: rotation 
curve 



1. Introduction 

The dynamical mass distribution in the Milky Way 
Galaxy is principally determined by analyzing the rota- 
tion curve on the assumption of circular rotation of the 
galactic disk (e.g., Sofue and Rubin 2001). In our previous 
work, a large scale rotation curve, called pseudo-rotation 
curve, has been obtained by combining the current rota- 
tion curve and radial velocities of globular clusters and 
satellite galaxies, which was used to model the Galaxy's 
mass components (Sofue, et al. 2009, Paper I; Sofue 2009, 
Paper II). While fitting to the inner rotation curves was 
satisfactory, the outermost mass structure related to the 
dark halo was still crude because of the large scatter of 
the used kinematical data. It has been shown that the 
kinematics of satellite galaxies is crucial for determining 
the global dark halo models. 

In this paper we construct a new large-scale, running 
averaged rotation curve, which we call the grand rotation 
curve. We analyze it more quantitatively and accurately 
compared with the previous work using the least squares 
fitting method. We adopt the de Vaucouleurs (1958) pro- 
file for the bulge, exponential disk (Freeman 1970), and 
NFW (Navaro, Frenk and White 1996) density profile for 
the dark halo. 

Throughout the paper, we adopt a galacto-centric dis- 
tance and circular velocity of the Sun of (i?o, Vb) = (8.0 kpc, 
200 km s" 1 ). The result of the local dark matter density 
will be compared with the most recent values from ex- 
tensive analyses of the galactic rotation curve (Xue et al. 



2008; Weber and de Boer 2010; Salucci et al. 2010). The 
fitted parameters for the bulge, disk, and dark halo will 
be used to estimate the baryon fraction inside the sup- 
posed Galaxy boundary at ~ 385 kpc, a half distance to 
M31. The estimated emission measure of hot gas will be 
compared with the observed X-ray background. 

2. Observed Rotation Curve and Velocity 
Dispersion 

Figure 1 shows the observed rotation velocities V(R) 
and galacto-centric radial velocities Vociobs, as repro- 
duced from Papers I and II, where we compiled kinemat- 
ical data from the literature and recalculated them for 
the galactic constants of (Ro,Vq) = (8 kpc, 200 km s _1 ). 
The inner data within R ~ 10 kpc are mostly from the 
traditional circular velocity measurements of disk objects 
near the Galactic plane, while some data with large scat- 
ter are the galacto-centric velocities of globular clusters. 
The data from 10 to ~25 kpc are the mixture of usual 
circular velocities and galacto-centric radial velocities of 
globular clusters. Those beyond R ~ 30 kpc are Vc&obs of 
distant globular clusters, satellite galaxies, dwarf galaxies, 
and member galaxies of the Local Group. See Paper II for 
detail and data references. We included the circular ve- 
locities from SDSS blue stars analysis by Xue et al. (2008; 
their model 1) by correcting for the systematic difference 
in velocities of about -20 km s _1 due to the our Vb = 200 
km s _1 instead of their 220 km s _1 . We also added the 
most accurate rotation velocities from VERA for stellar 
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Figure-included pdf ms is at URL http://www.ioa. s. u-tokyo.ac.jp/~sofue/htdocs/2012DarkHalo/ 



Fig. 1. [Top] Rotation velocities and galacto-centric radial velocities made from the data list in Paper II. Circles are disk objects, 
triangles are non-disk objects showing their galacto-centric radial velocities (not multiplied by \/2), and reverse triangles are circular 
velocities from Xue et al. (2008). [Bottom] Same, but in logarithmic scaling up to 1 Mpc. 



maser sources (Honma et al. 2007; Oh et al. 2009). 

A rotation curve is denned by circular velocity V cnc in 
a balance between the centrifugal force and gravity which 
is equivalent to the Virial theorem of a circularly rotating 
disk objects. For non-disk objects like globular clusters 
and satellites that are moving in statistical and/or random 
orbits in Virial equilibrium, we define a pseudo-rotation 
velocity as galacto-centric velocity Vqc corrected for the 
freedom of motion of individual objects: it is related to 
the observed galacto-centric velocity as Vqc = V% Vc&obs, 
where V% is the correction factor for possible tangential 
motion of the object as expected from its observed radial 
velocity. We thus define the rotation curve as 



W(R) = V ci 



(1) 



for disk objects (HI gas, CO gas, molecular clouds, HII 
regions, OB stars, and red giants in the disk), and 

W(R) = Vqc = V^Vbc:obB (2) 

for non-disk objects (globular clusters, satellites, dwarf 
galaxies, Local Group galaxies). 

In order to obtain a smooth rotation curve, we apply 
a running mean of the observed values in figure 1. We 
define a rotation velocity at Ri as 



V(IU 



^J_ N W(R 3 ) 



(3) 



27V +1 

where 2N + 1 is the number of used data points around 
R = Ri, which is taken every N steps in the data list in 
the order of increasing radius. In the present paper we 
take N = 5, so that the data arc running averaged every 
N = 5 points using their neighboring 2N + 1 = 11 objects. 
We also calculate the velocity dispersion by 



<r(Ri) 



W{Ri)] 



(2AT + 1)-1 



(4) 



Figure 2 shows the obtained rotation curve and ?? 
. Note that the observational errors of individual data 
points are much less than the scatter in this plot, 
i.e. cr(Ri) is much larger than the observational er- 
ror. In this meaning a(Ri) represents the real veloc- 
ity dispersion of the objects in each averaging bin. A 
table of the obtained rotation velocities used for fig- 
ure 2 is available from our URL: http://www.ioa. s.u- 
tokyo.ac.jp/~sofue/htdocs/2012DarkHalo/ 

3. The Models 

We assume that the rotation curve is composed of a 
bulge, disk, and dark halo components, and is calculated 
by 



V(R) 2 = V h {Rf + V d (R) 2 + V h (R? 



(5) 



where V(R), Vb(R), Vd(R), and Vh(R) are the circular 
rotation velocity at galacto-centric distance R, the corre- 
sponding circular velocity due to the bulge, disk and dark 
halo, respectively. Using the newly constructed rotation 
curve in figure 2, we construct a model rotation curve us- 
ing the least square fitting. 

3.1. Bulge 

The bulge is assumed to have the de Vaucouleurs (1958) 
profile for the surface mass density as 

1/4 
- I - 1 



E b (r) = E be exp 




(6) 



where k = 7.6695, £t>c is the surface mass density at the 
scale radius R = a b . Note that the scale radius a b is de- 
fined so that the enclosed mass (luminosity) within the 
cylinder of radius a b is equal to a half of the total mass 
(luminosity). The total mass is calculated by 



M b = 2tt 



-S b (r)c?r = rja 2 T, 



be • 



(7) 



with 77 = 22.665 being a dimcnsionless constant (Paper I 
for detail). The mass within a sphere of radius R is given 

by 

M b (i?)=4W p h (r)r 2 dr, (8) 



(9) 



where the volume density is calculated using E b as 



Ph{r) = - 



cffib(x) 



1 



dx 



-.dx, 



The circular rotation velocity is then given by 
V h (R)-- 



GM h 



3S(X), 



(10) 



where X = R/a^, and 
&{X) = 



Vx~ 2 



dxdy 



1/2 



.(11) 



In the fitting procedure, Mb and ab are taken as the two 
free parameters. The time consuming integral in the above 
equation was computed by calling a function subroutine 
which refers to a prepared table of SS[X). Figure 3 shows 
the 8$ function as compared with the other useful laws 
corresponding to unity total mass and unity scale length. 
Obviously, all laws tend to the Kclpcrial law at sufficiently 
large radii. 

3.2. Disk 

The galactic disk is approximated by an exponential 
disk, whose surface mass density is expressed as 
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Figure-included pdf ms is at URL http://www.ioa.s. u-tokyo.ac.jp/~sofue/htdocs/2012DarkHalo/ 



Fig. 2. Rotation curve V {Ri) of the Galaxy obtained by running mean of 11 = 2 X 5 + 1 neighboring data at every 5 data points in 
the order of increasing radius. The bars indicate least-square fit dispersion in each bin. Top panel shows the data within 25 kpc, 
and the bottom within the Local Group. Note the abscissa in the bottom figure is logarithmic. 



Figure-included pdf ms is at URL 
http:/ /www.ioa.s.u- 
tokyo.ac.jp/ ^sofue /htdocs / 2012DarkHalo / 



Figure-included pdf ms is at URL 
http:/ /www.ioa.s.u- 
tokyo.ac.jp/ ^sofue /htdocs / 2012DarkHalo / 



Fig. 3. The SS(X) function plotted against X = R/a h . For 
comparison, we show Keplerian law, Plummer potential, and 
flat disk 9{X) function with X = R/a d . 



R 

E d (i?) = E exp ( 



(12) 



where £o is the central value and ad is the scale radius. 
The total mass of the exponential disk is given by 



Ma = 



/>OG 

Jo 



(13) 



The rotation curve for a thin exponential disk is expressed 

by 



V d (R) 



GM A 



rid 



(14) 



where X = R/ad, and 
@(X) = (X/V2) x 

x [{/„ (X/2) K (X/2) - h (X/2) K t (X/2)}} 1/2 (15) 

with Ii and Ki being the modified Bessel functions 
(Freeman 1970). As the two free parameters, we chose 
the total mass Md and ad- 

3.3. Dark Halo 

For the dark halo, three mass models have been so 
far proposed: the semi-isothermal spherical distribution 
(Begeman et al. 1991), NFW (Navarro, Frenk and White 
1996), and Burkert (1996) models. In Paper II, we have 
shown that the isothermal model may be not a good ap- 
proximation to the observed rotation curve, and showed 
that the NFW and Burkert models may better fit the ob- 
servations. 

Since the NFW and Burkert profiles are essentially the 
same, except for the very inner region where the contribu- 
tion from the dark halo component is negligible, we here 
adopt the NFW profile. It is expressed as 
Po 



p{R) 



(16) 



x(i+xy 

where X = R/h, and po and h are the representative 
(scale) density and radius (core radius) of the dark halo, 
respectively. 

The enclosed mass within a radius 7? is given by 



M h (R) = inpoh* ( ln(l + X) - 



X 



l + X 



(17) 



The mass enclosed within a radius h is given by M£ = 
2A272poh 3 for X = 1, which we call the scaling mass of 
the dark halo. 



Fig. 4. [Top] Dispersion £ of the fitted result as a function of 
each parameter a^, a^, and h around its best value for other 
parameters being fixed to the best values. [Bottom] Same as 
figure 4, but for the masses M^, and Afj,. 

The circular rotation velocity at a given distance R is 
given by 



Vh(R) 



GM h {R) 



R 



(18) 



In the fitting procedure, we chose po and h as the two free 
parameters. 

4. Fitting Procedure 

We search for the best-fit parameters for the bulge, disk 
and dark halo using the least-squares method. The free 
parameters are Mb, Ob, -Md, ad, po and h. 

4-1. Data 

We divide the grand rotation curve and the velocities 
shown in table 2 and figure 2 into three parts: inner part 
for R < 5 kpc which is mainly used for the bulge fitting, 
middle part at 1 < R < i?2 = 40 kpc for the disk, and the 
outermost part at 1 < R < 400 kpc for dark halo fitting. 
In order to fit the rapidly varying bulge component, we 
added the innermost data at R < 1 kpc from the original 
data. 

The outermost boundary for the halo fitting was so cho- 
sen that it represents the gravitational boundary between 
the Galaxy and M31, or about a half distance, 385 kpc, 
to M31 at 770 kpc, as discussed in Paper II. 

4-2. The Least Squares Method 

We find the best fit values for Mj and ai for the bulge 
and disk by the least square method using the inner two 
parts, where the dark halo parameters are given as pro- 
visional initial values. Using the fitted bulge and disk 
parameters, we then find the best fit values of po and h 
for the dark halo using the outermost data. The thus 
found po and h are again used in finding the disk and 
bulge parameters, which is further used to find the halo 
parameters. This procedure is repeated several times until 
the fitted results reach sufficiently stable values. Below, 
we describe the procedure in more details. 

We define £ as 



Z[V ohs (Ri)-V(Ri)] 2 



(19) 



K-l 

where V a h& is the observed value plotted in figures 5 and 6, 
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Table 1. Grand Rotation Curve of the Galaxy^ Table 1. (continued) 



Radius 


Velocity 


Vclo. Dispersion 


Radius 


Velocity 


Vclo. Disp 


(kpc) 


(km s _1 ) 


(km s^ 1 ) 


(kpc) 


(km s _1 ) 


(km s _1 ) 


0.29 


241.68 


17.15 


6.05 


214.04 


5.47 


0.52 


242.67 


11.06 


6.13 


217.10 


11.67 


0.82 


226.67 


12.34 


6.21 


214.21 


15.33 


1.10 


214.93 


6.44 


6.28 


214.19 


18.46 


1.30 


210.32 


7.12 


6.35 


215.27 


13.26 


1.51 


203.26 


6.14 


6.42 


209.18 


21.97 


1.74 


201.01 


3.76 


6.51 


205.94 


20.53 


1.97 


199.71 


5.23 


6.61 


203.96 


18.62 


2.17 


196.65 


7.07 


6.71 


209.96 


1.64 


2.32 


193.97 


6.63 


6.79 


209.73 


3.30 


2.44 


192.61 


6.22 


6.87 


208.78 


3.44 


2.57 


191.40 


4.05 


6.96 


207.64 


2.83 


2.70 


191.84 


3.76 


7.06 


208.47 


4.98 


2.82 


192.93 


3.29 


7.14 


207.87 


4.96 


2.94 


192.51 


3.70 


7.21 


205.69 


2.87 


3.06 


192.54 


4.76 


7.28 


204.68 


2.82 


3.19 


197.59 


8.86 


7.35 


205.02 


3.67 


3.30 


200.61 


9.73 


7.41 


205.48 


3.38 


3.42 


200.16 


6.08 


7.48 


202.02 


7.30 


3.54 


199.11 


3.43 


7.54 


201.82 


7.20 


3.65 


200.92 


7.10 


7.60 


203.75 


2.28 


3.77 


204.13 


9.23 


7.66 


203.36 


3.02 


3.89 


208.48 


9.79 


7.72 


203.73 


3.58 


4.01 


215.06 


11.90 


7.78 


204.03 


3.10 


4.13 


217.24 


8.67 


7.83 


208.38 


10.80 


4.25 


212.75 


8.18 


7.87 


210.16 


10.21 


4.36 


210.44 


4.87 


7.92 


208.39 


7.09 


4.47 


211.66 


6.42 


7.95 


209.52 


11.52 


4.58 


207.17 


6.82 


7.98 


206.54 


10.40 


4.70 


204.11 


4.28 


7.99 


204.74 


4.87 


4.82 


206.15 


4.52 


8.04 


202.92 


7.16 


4.93 


209.06 


3.71 


8.17 


198.71 


16.40 


5.04 


209.52 


4.25 


8.38 


195.51 


18.80 


5.15 


207.54 


5.06 


8.61 


188.61 


13.70 


5.24 


208.90 


9.63 


8.88 


185.18 


13.73 


5.32 


208.38 


10.02 


9.16 


179.90 


15.12 


5.41 


205.36 


4.74 


9.35 


176.79 


47.92 


5.52 


205.82 


4.40 


9.43 


179.48 


49.56 


5.63 


210.14 


7.93 


9.51 


177.00 


19.42 


5.74 


210.94 


7.41 


9.61 


188.14 


48.17 


5.84 


208.27 


2.84 


9.84 


207.90 


59.99 


5.95 


210.60 


3.97 









f The data are running-averaged values using the com- 
piled data in Paper I and II (Sofue et al. 2009; Sofue 
2009), data from VERA (Honma et al 2007; Oh et al 
2009), and those from Xue ct al. (2008). Ascii data for 
this table is available from URL 

http:/ /www. ioa.s. u-tokyo.ac.jp/~sofue/htdocs/2012DarkHalo/ 
Note that, in order to avoid unnecessary smoothing of the 
steeply varying curve around the bulge component, the 
fitting was obtained by replacing the data in the table at 
R < 2 kpc with the denser original data points (before 
running average). 
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Table 1. (continued) 



Radius 
(kpc) 


Velocity 
(km s -1 ) 


Vclo. Disp 

(km s -1 ) 


10.13 


203.08 


61.44 


10.33 


189.00 


53.36 


10.43 


177.57 


30.87 


10.56 


176.56 


26.96 


10.81 


185.25 


72.22 


11.09 


198.74 


70.78 


11.31 


190.52 


31.04 


11.59 


209.52 


64.58 


11.89 


218.16 


69.56 


12.18 


185.67 


60.42 


12.46 


177.61 


53.11 


12.76 


187.25 


41.79 


13.16 


183.61 


63.00 


13.63 


169.95 


60.50 


14.33 


183.58 


58.46 


15.35 


210.68 


64.01 


16.57 


216.70 


49.18 


18.43 


175.61 


69.31 


20.43 


181.42 


91.02 


23.81 


199.71 


92.24 


30.42 


164.53 


67.49 


42.97 


134.85 


66.69 


61.88 


142.41 


44.90 


81.39 


145.05 


49.65 


125.49 


120.67 


82.77 


353.29 


82.83 


83.95 


647.55 


105.77 


96.74 


817.71 


118.92 


80.17 


1012.54 


126.69 


98.22 


1279.93 


140.53 


102.43 



and K is the number of used data points in the individual 
data parts. Therefore, K is the number of plotted points 
at R = to Ri for the fitting of the bulge, from R\ to i?2 
for disk, and from i? 2 to Ry lm for the dark halo. 

We here regard all the plotted points to have the same 
weight for the following reason: Since the plotted values 
in figure 2 are the running means of the observations, the 
bars in the plot are not measurement errors. Particularly, 
the outer data points include real velocity dispersion, 
while we are interested in the mean values. In order to 
discuss the dark halo using the outermost velocity infor- 
mation, it is not practical to apply such weighting that is 
proportional to the inverse of the square of measurement 
error bars, which should give a very small weighting to 
the outermost rotation curve. In order avoid unnecessary 
smoothing of the steeply varying curve around the bulge 
component, the fitting was obtained by replacing the data 
in the table at R < 1 kpc with denser original data points 
before running average. 

Since the number of parameters to be determined is six 
(ai, Mi, 02, M2, po, and h), it is not practical to obtain 
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the best fit by one time iteration, we apply a step-by-step 
iteration as the following. 

We first use the inner rotation curve, and find param- 
eters that represent the bulge and disk components. The 
dark halo density is set to an arbitrary value. Thereby, we 
find values of a\ and M\ that minimizes £ for an appropri- 
ately assumed initial values of a-i and M 2 . We calculate 
£ by changing a\ and M\ at small steps, and fix the best 
fit values that give minimum £. Thus found a\ and M\ is 
used to find the best fit 02 and M2. This procedure is re- 
peated until the best fit values both for the bulge and disk 
reach stable values. The fixed values of 01, 0,2, Mi and M2 
are used to find the best fit values of h and po by the same 
procedure using the outer rotation curve at 1 < R < R\i m 
kpc. Using the best-fit halo values, the bulge and disk fit 
is repeated. By repeating these procedures several times, 
we obatin a stable set of the six parameters. 

4-3. Variations of the fitted parameters 

Figure 4 show the variation of the calculated dispersion 
£ around the best fit value of scale radius at>, ad, and h as 
a function of each parameter with the other parameters 
fixed to the best values. Figure 4 shows the same but for 
the mass parameters. For the dark halo, the abscissa is 
taken as the scaling dark mass M^ . 

4-4- Error estimates 

The errors of individual fitted parameters are evaluated 
as the range of the parameter that allows for an increase 
of the dispersion £ by a factor of 1.1, i.e. 10% increase of 
the dispersion. Errors in quantities which are calculated 
using the fitted parameters are computed by considering 
propagation through the corresponding equations. 

5. The Best-Fit Galactic Parameters and 
Rotation Curve 

Table 5.2 shows the best-fit parameters for individual 
mass components of the Galaxy, and figures 5 and 6 show 
the calculated rotation curves from these parameters. 

5.1. Radii and Masses of the Bulge and Disk 

The scale radius, i.e. the half-pro jected-mass radius, 
a h = 522 pc, and the total mass M b = 1.652 x 1O 1O M of 
the bulge are consistent with those given in the literature. 
The disk scale radius ad = 3.19 kpc is within the range of 
often quoted values, and the disk mass of is Md = 3.41 x 
10 10 M o . The bulgc-to-disk mass ratio is M b /M d = 0.48. 
This ratio is closer to that for Sbc galaxies: it is 1.0 for Sb 
and 0.38 for Sbc galaxies (Koeppen and Arimoto 1997). 

5.2. Dark halo radius and density 

The best fit scale radius of the dark halo was obtained 
to be 

h= 12.53 ±0.88 kpc, (20) 
and the scale (representative) density 

po = (1.06 ±0.14) x 10- 2 M Q pc- 3 (21) 
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G 
or 

po = 0.403 ±0.051 GeV cm" 3 , (22) 

which give the least £ value of 

£ = 9.6kms _1 (23) 

The local value near the Sun at R = 8 kpc is estimated 
to be 

pf = (6.12 ± 0.80) x 10" 3 A/ Q pc" 3 , (24) 

or 

p% = 0.235 ±0.030 GcV cm" 3 . (25) 

6. Discussion 

6.1. The Grand Rotation Curve 

We have constructed the grand rotation curve of the 
Milky Way Galaxy using the data compiled in Paper I 
as well as by adding the most recent accurate measuring 
results. The grand rotation curve cover covers a wide area 
from the Galactic Center to the Local Group space. The 
curve well coincides with the circular velocity curve by 
Xue et al. (2008) from 7.5 to 55 kpc. 

By the least-squares fitting to the thus obtained ro- 
tation curve within the boundary of the Milky Way at 
R ~ 400 kpc, the parameters of the three mass compo- 
nents, the de Vaucouleurs bulge, exponential disk and 
dark halo with NFW profile, were determined at high ac- 
curacy. 

The fitting was obtained more quantitatively and ac- 
curately compared to our earlier works (Paper I, II) by 
removing the artificial assumption of the critical radius at 
which the dark halo contribution becomes equal to that 
from the disk. The resulting fitting parameters are listed 
in table 5.2, and the fitted rotation curve is shown in fig- 
ures 5 and 6. 

6.2. The Best-Fit Rotation Curve 

In the linear representation in figure 5, the calculated 
rotation curve well reproduces the observed data, except 
for the wavy and bumpy behaviors, which are attributable 
to local rings and/or arm structures as discussed in Paper 
I. The fitting seems satisfactory almost equally both in the 
two cases with and without including the dark halo. In 
other words, it is difficult to discuss the dark halo problem 
using the current rotation curves up to R ~ 20 — 30 kpc. 

Obviously the logarithmic expression in figure 6 makes 
it much easier to discriminate the dark halo contribution. 
Comparison of figures 5 and 6 demonstrates how the grand 
rotation curve is efficient and essential to discuss the dark 
halo in the Milky Way. The NFW model appears to re- 
produce the observations sufficiently well within the data 
scatter. 

Another interesting fact is that the most inner part at 
R < 0.2 kpc systematically deviates from the observations: 
the observations more rapidly decrease toward the nucleus 
than the model. This indicates that the de Vaucouleurs 
law may not be the best expression of the innermost mass 
distributionc. 
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6.3. The Galactic Parameters 

Table 5.2 summarizes the obtained parameters by the 
present analysis, and figures 5 and 6 show the fitted ro- 
tation curves. The parameters for the bulge is not much 
different from the current works, because the correspond- 
ing rotation profile at R ~ 1 kpc has a steep peak and the 
least square fitting is more independent and effective com- 
pared to the fit to other extended components like the disk 
and halo. The total mass of the bulge and disk system, 
Mb+d = M h +M d = 5.06 x 10 10 M Q , is larger than the dark 
mass M h (R < R Q ) ~ 2.75 x 10 10 M o inside the solar circle. 
The bulge-to-disk mass ratio was found to be 0.48 is closer 
to that for Sbc (0.38: Koeppen and Arimoto (1997)). Our 
Galaxy may be, thus, a slightly late Sb galaxy. 

We here confirm that the total mass of the bulge, disk 
and the dark halo within the solar circle Afb+d + M^{R < 
Rq) = 7.8 x 10 10 Mq ~ Mq. is comparable to the simple 
spherical mass estimate for Vq — 200 km s _1 at Rq = 8 
kpc, Mq = R V 2 /G = 7.44 x 10 10 A/ Q . The presently larger 
value is mainly due to the disk effect, and partly due to 
extended bulge and disk components beyond the solar cir- 
cle. 

Xue et al. (2008) obtained total mass of 4.0 ± 0.7 x 
10 11 M Q within R = 60 kpc, which corresponds to 3.6 x 
1O U M for V Q = 200 km s" 1 . In our fitting result, the total 
mass within 60 kpc is calculated to be 2.4±0.4 x 1O 11 M , 
and is consistent with their value within the error. 

6.4- NFW Dark Halo and Local Dark Matter Density 

The NFW profile was found to fit the grand rotation 
curve quite well. The declining part in the outermost ro- 
tation velocity predicted by the NFW profile was clearly 
detected for the first time, showing the good fit to the out- 
ermost grand rotation curve at R ~ 40 — 400 kpc covering 
the Galaxy's gravitational boundary. 

The local value of the dark matter density in the Solar 
neighborhood as derived from the present analysis, p® = 
0.235 ±0.030 GeV cm" 3 , is slightly smaller than the re- 
cently published values 0.43 ±0.10 GeV cm _3 by Salucci 
et al. (2010), while consistent with the value0.2 - 0.4 GeV 
cm _3 by Weber and de Boer (2010). 

6.5. Baryon Fraction 

The total Galaxy mass enclosed in the supposed grav- 
itational boundary of the Galaxy at 385 kpc, a half dis- 
tance to M31, is ~ 7.03 x 10 11 M Q . As discussed in Paper 
II, this mass is not sufficient to gravitationally stabilize 
the Galaxy and M31 system, and therefore the two galax- 
ies may be embedded in a larger-scale dark matter system 
properly constructing the Local Group. 

If we compare the total dark matter mass with the bulge 
and disk mass, we obtain an upper value to the baryon- 
to-dark matter ratio in the entire Galaxy shared by the 
bulge and disk to be 

/bd = M b+d /(M b+d + M h ) ~ 0.072 ± 0.018. (26) 

Athough this fraction is much smaller than the WMAP's 
cosmic value of baryon fraction, 0.17 (Dunklcy et al. 
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No. 



Dark Matter Halo of the Milky Way Galaxy 

Table 2. Best-fit parameters for the mass components of the Galaxy 



Mass component 



Mass; Density 



Scale Radius 



Dispersion £ 



Bulge parameters 
Disk parameters 
Bulge+Disk Mass 
Bulge/Disk Mass ratio 



M b = (1.652 ±0.083) x 10 10 M Q 
M d = (3.41 ± 0.41) x 1O 1O M 
M h+d = (5.06 ±0.97) x 10 10 M Q 
M b /M d = 0.48 ±0.09 



a h = 0.522 ±0.037 kpc 
a d = 3.19 ±0.35 kpc 



8.63 km s" 1 
9.53 km s" 1 



Dark halo parameters 

Local DM densityt 
at Sun (R = 8 kpc) 


po = (1.06 ±0.14) x 10- 2 Af o pc- 3 
= 0.403 ±0.051 GeV cm" 3 

pf = (6.12 ±0.80) x 10- 3 M Q pc- 3 
= 0.235 ±0.030 GeV cm- 3 


h= 12.53 ±0.88 kpc 


9.57 km s" 1 


Dark halo mass^ 


M h (R < 8kpc) = (2.71 ± 0.42) x 1O 1O M 
M£(< h) = (5.05 ±0.78) x 1O 1O M 
M h (< 20kpc) = (8.87 ±1.37) x 1O 1O M 
M h = M h (< 385kpc) = (6.52 ± 1.01) x 1O U M 


(Scaling DH mass) 




Galaxy Mass 


Mb+d+h = (7.03 ± 1.01) x 1O U M 


(R < 385 kpc) 




Baryonic fraction 


M h+d /(M h+d+h ) = 0.072 ± 0.018 







f 1 M pc 3 = 38.16GeVcm 3 . $ Error propagation included. 

Figure-included pdf ms is at URL http://www.ioa.s.u-tokyo.ac.jp/~'Sofue/htdocs/2012DarkHalo/ 

Fig. 5. Least-squares fit by bulge, disk and dark halo to the grand rotation curve. Thick line represents the fitted rotation curve, 
and thin lines show individual contributions from bulge, disk and halo. Observed velocities V(Ri) are shown by open circles, and 
the most recent accurate results from VERA (Honma et al. 2007; Oh et al. 2010) are shown by squares. 

Figure-included pdf ms is at URL http://www.ioa.s. u-tokyo.ac.jp/~sofue/htdocs/2012DarkHalo/ 

Fig. 6. Same as figure 5, but in logarithmic scaling displaying up to 1 Mpc. The dark halo is satisfactorily fitted by the NFW 
model. This figure demonstrates that a grand rotation curve up to ~ 1 Mpc is essential in order to quantitatively analyze the dark 
halo. 



2009) , it may be compared with the mean baryon frac- 
tion, ~ 0.12, observed for groups of galaxies (Andreon 

2010) . If we assume that the fraction in the Galactic 
halo of radius of ~ 385 kpc is the same, a small fraction, 
~ 0.05 (= 0.12 — 0.072), corresponding to baryon mass of 
~3x 1O 1O M , may possibly exist in the form of gas filling 
the dark halo at density of ~ 6 x 10~ 6 cm -3 . 

If the gas is diffuse HI or H2 gas, the column density 
would be ./Vh ~ 7 x 10 18 H cm -2 and the velocity width 
as wide as ~ 100— 200 km s _1 . The brightness tem- 
perature will be ~ 0.02 K for HI or - 0.2 mK for CO 
( J = 1 — 0). Such low brightness and diffuse line emis- 
sions are extremely difficult to be detected. Moreover, it 
is unlikely that such diffuse gas can remain neutral, be- 
cause the internal motion is highly supersonic, and the 
gas is easily ionized. Alternative forms may HI and/or H2 
clouds, but it is also unlikely that such a large amount of 
HI or molecular clouds, almost ~ 10 2 times more massive 



than the interstellar gas in the galactic disk, have not been 
detected by the current radio surveys (e.g. Wakker and 
van Woerden 1997) . 

One possible way for the gas to coexist in the dark 
halo would be that they are in the form of X-ray emit- 
ting hot plasma. The temperature will be ~ 10 6 K or 
~ 1 keV, and electron density n c ~ M bary / (mH47ri? 3 /3) ~ 
0.6 x 10 _5 cm -3 , corresponding to emission measure of 
~ x 10 _5 pc cm 2 for the depth of R ~ 385 kpc. This emis- 
sion measure is much smaller than the observed value 
of ~ x 10 _2 pc cm 2 toward the North Galactic Pole by 
ROSAT (Sidher et al. 1999; McCammon et al. 2002). 
Hence, it is possible that the apparently 'missing' baryons 
are in the form of hot halo gas, and are already 'seen' in 
the ROSAT X-ray background emission, sharing its small 
portion. 
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